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Introduction. Suppose that tp and ip are functions defined on the real line and let 



x-y 

where x is the characteristic function of some union of intervals 

J = Uf =1 [a 2j -i,a 2j ). 
We regard K(x,y) as the kernel of an integral operator K 



> 

m 

O 

O' Kf{x)= I K(x,y)f(y)dy 

Jr. 



For special choices of ip and ip, Fredholm determinants of the form 



gn! det(l-K), 



arise in the theory of random Hermitian matrices [8] . In this theory it is a result of 
Gaudin and Mehta that such determinants give the probability that no eigenvalues 
for an N x N random matrix lie in the set J (the limit iV — > oo is also of interest 



in the context of various scaling limits [12]). More generally the derivative 

(-l) n d n 



det(l - \K) 



A=l 



is the probability that exactly n eigenvalues of the random matrix lie in the interval 
J. 

Fredholm determinants of a similar sort arise for certain solvable models in two 
dimensional quantum field theory [4], [3]. In the pioneering work [4], it was shown 
for the simplest N —>■ oo limit in random matrix theory and for impenetrable Bosons 
in two dimensional field theory, that it is possible to express the relevant Fredholm 
determinant in terms of the solution to a completely integrable nondinear system of 
differential equations. In [3] some other field theory models are shown to have this 
same property. In reference to this phenomenon, Its, Izergin, Korepin, and Slavnov 
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refer to kernels of type (0.1) as "completely integrable integral kernels". Recently, 
Tracy and Widom [15], have given a systematic treatment in the context of random 
matrix models. They show that if the functions ip and ip satisfy a linear differential 
equation 

(0.2) 
where 

A(z) B(z) ' 
-C(z) -A(z)_ ' 

and m(z), A(z), B(z) and C(z) are polynomials, then the Fredholm determinant, 
det(l — K), can be expressed in terms of the solution to a system of non linear 
differential equations. 

In this note we will show that the differential equations found by Tracy and 
Widom are special cases of monodromy preserving deformation equations, and the 
analysis of the Fredholm determinant is a special case of the general analysis of 
r — functions [5]. This is not the first time that a connection between random 
matrices and monodromy preserving deformation theory has been observed. Moore, 
[9], [10], has already observed that families of differential equations of type (0.2) 
play a role in the analysis of the partition function for random matrix models. In 
his work the singularity at infinity for (0.2) is important. The deformations of the 
coefficients m, A, B and C that he considers all preserve the "Stokes' multipliers" for 
(0.2), i.e., the monodromy data at infinity. The partition function is a r — function 
for the resulting monodromy preserving deformation, and the general connection 
between r — functions and non-linear equations accounts, in particular, for the 
expression which gives the partition function for two dimensional gravity in terms 
of a Painleve transcendent. Some interesting mathematical problems that arise in 
Moore's analysis have been solved by Beals and Sattinger in [1]. In particular they 
prove that certain deformations of (0.2), governed by compatible linear differential 
equations for <p and ip, do indeed leave the Stokes' multipliers at infinity invariant. 

Here a local perturbation of (0.2) with simple poles at the points aj will play a 
central role (see (1.14) below). The compatibility of the linear differential equations 
(1.14) and (1.15) will lead naturally to the non-linear differential equations found by 
Tracy and Widom. Realizing the non-linear equations as compatibility conditions 
for linear systems is not new. In [2] Harnad, Tracy and Widom obtain the differential 
equations for the sine kernel and the Airy kernel as compatibility conditions for 
linear systems that are essentially the same as the linear systems we find below. 
Apparently, these linear systems were guessed on the basis of experience with the 
dynamical systems on coadjoint orbits of loop algebras that are the central theme 
of [2]. One advantage of our treatment is that the linear systems are very simple 



d 

dz 



= Q(z) 



Q(z) = —- 
m(z) 
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to get at and the connection with dynamical systems on coadjoint orbits of loop 
algebras is then direct (in the setting which is explained in [2]). 

This paper is organized in two parts. In the first part we show how to realize the 
Tracy- Widom equations as monodromy preserving deformation equations, staying 
close to the integral equation which is central in their work. In the second part 
we introduce the formalism of Malgrange [7] in order to discuss the deformation 
theoretic r — function. Up to a trivial factor we show that this is the determinant, 
det(l — XK). The formalism of Malgrange is, incidently, better adapted to deal with 
analyticity questions in the aj parameters. 

The author would like to thank C. Tracy and H. Widom for helpful conversa- 
tions about their work. 
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§1 The Deformation Problem 

We will now introduce the deformation problem which will lead to the linear equa- 
tions (1.14) and (1.15). We will consider the problem in light of the formalism 
for Cauchy-Riemann operators introduced in [11]. It is not essential to understand 
the formalism in [11] to follow what goes on but it does provide some motivation. 
Suppose to begin that 

(p(z) 

is a solution to the differential equation (0.2) defined for z G C, the complex plane. 
Here m, A, B, and C are assumed to be polynomials as mentioned above. In case m 
is a non constant polynomial it may happen that the differential equation (0.2) does 
not have any globally defined single valued solutions. To avoid this complication at 
the start we will suppose that m = 1. Let 



(p(z) 



be a second independent solution to the same differential equation (0.2). Since 
the trace of the matrix on the right hand side of (0.2) is the determinant of the 
fundamental solution 

ip ip 



is constant and we can choose p and ip so that this determinant is 1. 
consider the Cauchy Riemann operator 



We wish to 





-l 


~d z 


" 




> £" 


ip ip 







d z _ 




Ip Ip 



(1.1) 



In this formula d z = \{d x + id y ) as usual. Of course, since ip, cp, ip, and ip 
are analytic functions in C this C-R operator acts precisely as the standard C-R 
operator 

~d z 
d z 



(1.2) 



The difference is that we will take the domain of (1.2) to be the Sobolev space 
if 1 (C) so that (1.1) represents an operator with with domain 



ip <£ 
1p Ip 



-i -1 



^(C). 
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Note that since the determinant of the fundamental solution is 1 the inverse 



'if £~ 


-1 


lp -If 


lp 1p 




_-<0 tp _ 



Next we wish to introduce a different domain for the C-R operator (1.1) which 
incorporates functions with fixed branching behavior along an interval [ai , a 2 ] on 
the real axis. It will be obvious how to incorporate more than one interval [cij, a^+i] 
in the formalism so I will confine my attention to just one interval for simplicity. 
We wish to include in the domain of our new C-R operator functions which change 
by a factor of 

"1 2ni\~ 
1 



M = 



in a counterclockwise circuit of the left endpoint, ai, and by the inverse of this 
matrix in a counterclockwise circuit of the right endpoint, a 2 - To understand the 
resulting C-R operator we would like to realize it as a similarity transform 



d z 
d z 



Y 



of the standard C-R operator by a multivalued analytic function Y(z). So that the 
domain of this operator is the same as the domain of (1.1) near infinity we want 



Y(z) ~ F(z) for z near 00, 



(1.3) 



where we write 



F(z) 



ip(z) ip(z) 
ip(z) i/j{z) 



To find Y(z) we split the problem into two pieces. Let Y + (z) = Y(z) for z E {z : 
z 2 > 0} and Y_{z) = Y(z) for z E {z : z 2 < 0}. Then a more precise formulation of 
(1.3) is that 



FY1 



1 
1 



(1.4) 



where H± is the usual Hardy space for the upper (+) and lower (-) half planes. On 
the real axis the boundary values Y + and Y"_ are related by 

Y+ = Y_((I- X )+XM) 

where x = X[oi,a 2 ] is the characteristic function of the interval. In order that the 
integral operator K (rather than its transpose) arises in the fundamental integral 
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equation (1.6) below it will be convenient to temporarily work with Y ± . The 
monodromy condition relating the boundary values of Y^ 1 becomes 

Y- 1 = ((I-x) + xM)Y-\ 

We can rewrite this last equation as 

YZ 1 - F' 1 = Y' 1 - F' 1 + X (M - ^Y' 1 . (1.5) 

Let P + be the usual orthogonal projection on H + and define 

pF = F -lp +F 

Then Pf is the analogue of P + for the C-R operator (1.1). We can write (1.4) as 

Pf (YT 1 -F- 1 ) =0, 

P+ (Y+ 1 - F' 1 ) = Y' 1 - F~\ 

where we are not multiplying operators anymore but just applying Pf to the 
columns of the matrices Y± — F~ l . Applying Pf to the columns of (1.5) and 
making use of these last relations one finds 



r+ 1 + P$ X (M - I)Y^ = f- 1 . 



;i.6) 



Since the kernel of P + is — ^(x — y) 1 (with x = x + iO) we find that the kernel 
of the operator P^%(M — I) is 



x-y 

Now let K denote the operator with kernel 

x(y) 



ip(x)ip(y) - (p(x)%j){y) 
<p(x)ifj(y) - ip{x)ip{y) 



x-y 



(<f(x)4>{y) - 4>(x)(p(y)) . 



This is a non-singular kernel. Let K denote the operator with kernel 

U(x)<p(y) - (p(x)^y)) • 
x — y V / 



This is a singular kernel. Because of the connection with P + it is understood that 
the variable x should be allowed to approach the real axis from above (i.e., x + iO). 
In terms of K and K the equation (1.6) becomes 



1 -XK 
1 - XK 



y + 1 



(1.7) 



This is solvable provided 1 — XK is invertible. We suppose that this is the case in 
what follows. As we will see in the second section, the r — function for this problem 
is just det(l — XK). 

Solving equation (1.7) for Y^ 1 one finds 



1> 



- XKP 
-P 



-(p + XKQ 

Q 



where 



Q=(l- XK)- 1 ^, 
P =(1 - AiTpV- 

We have here adopted the notation from [15] to facilitate comparisons. From the 
integral equation 

I - P + (F X (M - ^Y- 1 ) , (1.8) 



FY' 1 



it is clear that FY + 1 



has an analytic continuation into the upper half plane. Since 
F is analytic and invertible in the whole plane, Y + will be analytic in the upper 
half plane. The same argument shows that Y_ has an analytic continuation into the 
lower half plane. The two functions Y+ and Y_ fit together to give a single analytic 
function Y(z) which has a branch cut joining a\ and a^- We want to analyse the 
singularities of Y(z) in a neighborhood of the points aj and near oo. First we 
consider the local singularity at a\. Define f(z) = log(z — ai), with the branch cut 
for the logarithm taken to be the ray [ai, oo). Apart from the location of the branch 
cut it does not matter which determination of the logarithm we use. We claim that 
the matrix 

1 



Y(z) 



G(z), 



(1.9) 



is single valued and analytic in a neighborhood of z = a\. The left hand side of 
(1.9) is single valued because the monodromy of the matrix 



A/(*) 
1 



as z makes a circuit of a± exactly cancels the monodromy that is built into Y(z) 
by construction. The left hand side of (1.9) will have a removable singularity at a\ 
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because P(z) and Q(z) are regular at z = a\ and K introduces only a logarithmic 
singularity at z = oi. 

Taking determinants of both sides of (1.9) shows that det(Y(z)) is an analytic 
function of z near 2 = ai. In a similar fashion one can show that det(Y(z)) is 
analytic near z = a 2 - From (1.8) it follows that 

det( J P" 1 y(z)) = det(Y(z)) = 1 + 0(z _1 ), 

for z in a neighborhood of 00. Thus det(Y(z)) = 1 for all z. It follows from this 
and the earlier expression for Y that 



Y 



Q <£-\KQ 
P ip- XKP 



which displays the analogous roles played by (p, and ip and Q, and P. 
Differentiating (1.9) with respect to z one finds 



XY(z) 



1 




z — CL\ 

which may be easily rewritten 



+ Y'(z) 



1 Xf(z) 
1 



= G'(z), 



Y'(z)Y-\z) = 



X 



z — d\ 



-Y{z) 



1 




Y-\z) + G'(z)G- 1 (z). 



(1.10) 



Now using the fact that the determinant of Y is 1 and the formula for Y in terms 
of P and Q above one finds that 



Y'(z)Y- 1 (z) 



X 



z — d\ 



Q1P1 -Ql 
Pi -Q1P1 



+ O(l) for z ~ Oi. 



(1.11) 



In this formula we've written qj = Q(aj) and pj = P(aj). In a precisely analogous 
fashion we calculate the residue of Y'Y~ X at z = a 2 to find 



Y'(z)Y- 1 (z) = 



X 



z - a 2 



Q2P2 

p\ 



-4 

-Q2P2 



+ 0(1) for z ~ a 2 . 



(1.12) 



Now define 

A=-P+F(M-I) X Y+ 1 . 
Differentiating (1.8) with respect to z one finds 

F'Y' 1 - FY^Y'Y' 1 = A', 
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or 



F'F~ 1 (FY~ 1 ) - (FY~ 1 )Y'Y~ 1 = A'. 
Since FY -1 = I + A this can be written 

Y'Y' 1 = (J + A)" 1 0(/ + A) - (/ + A) _1 A'. 

Since A(z) = 0(z~ 1 ) and A'(z) = 0(z~ 2 ) it follows that the principal part of 
FT" 1 at oo is the same as the principal part of 



(/ + A)- 1 0(/ + A). 



(1.13) 



It now follows from this observation and (1.11) and (1.12) that there is a matrix 
valued polynomial M(z) of the same degree as Q(z) so that 



^Y- 1 = M(z) + ^ 

dz Z — d\ z — 0,2 



(1.14) 



where 



QkPk ~q% 
v\ -QkPk 



Now we determine the differential equation satisfied by Y in the aj variables. 
Differentiating (1.8) in the a variables one finds 

d a YY- x = Oiz' 1 ) as z -> oo. 

Differentiating (1.9) in the a variables one finds the principal part of d a YY~ x at a,j 
is: 

(-l) k ^da k . 

z- a k 

Combining these two observations one finds that the function d a YY _1 is analytic 
in z on P 1 except for simple poles at a\ and a<i- It has a zero at z = oo and so 



daYY- 1 



XCi AC 2 

-da\ H da<2- 



z — d\ 



z - a 2 



(1.15) 



Now we apply these observations to deduce non-linear deformation equations. If 
one examines the equality of mixed partials 



da,j dz 



dz da^. 



using (1.14) and (1.15), one obtains an equation relating rational functions of z with 
poles at cij and oo. Equating the principal parts at the finite poles aj, one finds 



d a d = [M(a 1 ),C 1 ]da 1 + X[C U C 2 
d a C 2 = [M(a 2 ),C 2 ]da 2 + X[C u C 2 



d(a 2 - at) 
a 2 — a\ 

d(a 2 - at) 
a 2 — ai 



(1.16) 



where M(aj) is the value of the matrix polynomial M(z) at z = aj. Equating the 
principal parts at oo one finds 



dM(z) 
daj 



(-iyx 



M(z) - M(a 3 



z — a n 



-, Cj 



(1.17) 



Note that in this last equation each side is a polynomial function of z of order 1 
less than the degree, n, of M(z) (since the highest order term in M(z) is clearly 
independent of cij). In particular if we write 

M(z) = M n z n + M^z 71 ' 1 + ... 

then it is easy to see that (1.17) implies 



(-iy\[M n ,Cj 



In order to make (1.16) and (1.17) more explicit we must use (1.13) to compute the 
polynomial M(z). To make use of (1.13) we must in turn calculate the asymptotics 
of A(z) at z = oo. Substituting for Y + and F in the formula for A(z) one finds 



A(z) 



Ax(y) 



-<p(y)P(y) <p(y)Q(y) 
-^(y)P(y) ^(y)Q(y) 



Ir z-y 

Expanding the denominator in a geometric series 



dy. 



1 



1 



z-y z 



y 



=-!+-+ 



one finds that 



A(z) 



A 



-v 


u 


A 

+ — 


-vi 




_-w 


v _ 


z z 




^1 . 



+ 0(z" 3 ), 



(1.18) 
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where 



Vi — 



Vn = 



Wj 



x(y)y J v(y)Q(y)dy, 



R 



/ x(y)y 3 v{y)P{y)dy, 

Jr 

x(y)y J i>(y)Q(y)dy, 
x{y)y J ^(y)P(y) dy, 



(1.19) 



R 



and the notation is again taken from [15]. We use (1.18) and (1.13) to calculate the 
matrix M{z) in a few of the cases considered in [15]. 
Our first example is the sine kernel with 

(f{z) = sin(z), 
ip(z) = cos(z). 



In this case 

n{z) = 

and a simple calculation shows that 



1 
-1 



M SineW = 



1 
-1 



Our second example is the Airy kernel with 

(p(z) = Ai(z), 
iP(z) = Ai'(z). 



In this case 



Q(z) 



1 

z 



Without difficulty one may use (1.13) and (1.18) to compute 

M Mrj (z) 



-Xu 1 
z — X(v + v) Xu 



(1.20) 



(1.21) 



where we've written u = uq, v = vq and v = vq. For A = 1 this simplifies somewhat 
since vq = vq in this event. 
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Our last example is one that arises in 2D quantum gravity matrix models 
(the case of pure gravity). A function is given which satisfies the differential 
equation 

?"(x) + 6Z{x)e(x)+4 = 0, 

which is a special case of the string equation. The functions (p(z,x) and ip(z,x) 
now depend on the parameter x and satisfy the differential equation 



d_ 

dz 



= n{z) 



where 



n(z) = 



z+^(x) 



- z *-u{ X ) Z -\e{x)- l A"{x) \?{x) 



This is incidently, an example of a family of differential equations in which the 
Stokes' multipliers at oo are independent of x [9], [10]. Again one may use (1.13) 
and (1.18) to compute 



Mi 



Gravity(*) = fi (*)+ M i* + M o 



where 



M x = 



Xu 1 

X(v + v) —Xu 



and 



M = A 



— Xuv + 7}£(x)u — w + Ui 



Xu 2 + (v + v) 



~S i 

A (uw — v(v + v)) + £(x)v + Xuv — -^^(x)u + w — u\ 



where we've written u = uo, v = vq and v = vq. 

It is straightforward to use these results to determine the differential equations 
(1.16) and (1.17). We will not do this here but instead we will refer the reader to 
[15], where the non-linear equations for qj, pj, uj, ... are written out, and some 
interesting examples beyond those mentioned above are considered. Tracy and 
Widom also show that equations (1.16) and (1.17) sometimes have additional first 
integrals which can be used to reduce the number of variables in the deformation 
equations. In several examples they show that the the deformation equations can 
be integrated in terms of Painleve transcendents (thus generalizing the original 
observation in [4] that Painleve V transcendents arise in the integration of the 
deformation equations associated with the sine kernel). 

To conclude this section we will describe the changes which need to be made 
when the differential equation for <p and ip already has poles in the finite plane. The 
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simplest example of this sort considered by Tracy and Widom is the Bessel kernel, 
for which 

<p(z) = J a (y/z), 
^{z) = zJ' a {V~z). 

One can see from this that ip(z) and ip(z) will not be single valued in the complex 
plane. One can deal with this situation by introducing a branch cut for tp and ip 
joining to oo. Since the intervals [a,j, a,j+i] in this model will all be subsets of 
the positive real line it is convenient, in this case, to choose the branch cut on the 
negative real axis. The differential equation satisfied by <p and ip is 



d 




1 


1" 






dz 


A 


z 









which has a simple pole at in addition to an irregular singular point at infinity. 



Extending 



to a matrix valued fundamental solution F with determinant 1 is 



possible as before provided one admits a branch cut for F on the negative real axis. 
Let F + denote the analytic function F in the upper half plane and let F_ denote 
the analytic function F in the lower half plane. The boundary values of these two 
functions on the real axis are related by 



F+ = F_ (I + xo(M - I)) 



(1.22) 



where xo is the characteristic function of the interval (— oo,0] and M is the mon- 
odromy matrix associated with the fundamental solution F. We will want the 
monodromy of the fundamental solution Y to match that of F on the negative real 
real axis and so the relation we desire between the boundary values Y± is 



Y ~ =(! + Xo(M -I) + X (M - I)) Y~ 



-l 



(1.23) 



Using (1.22) to write F_ in terms of F + , and the fact that xoX = 0, one can multiply 
(1.23) on the left by F_ to get 

F_YZ l -1 = F + Y~ l -1+ F +X (M - I)Y~ 1 . 

If we now insist that 



P± (F±Y± X -I)= F±Y± 1 - /, 



as before, then we find 



F+Y+ L = I - P + F + (M - l) X Y+ L , 



(1.24) 
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instead of (1.8). Now let 



A = -P+F(M-I) X Y- 1 

and differentiate (1.24) with respect to z to obtain 

1 = (/ + A) _1 0(/ + A) - (/ + A) _1 A'. 

Now since A(z) = 0(z~ 1 ) and A'(z) = 0(z~ 2 ) we can conclude that the principal 
part of zY+(z) Y + (z) _1 at z = oo is the same as the principal part of 



{I + A(z))- 1 Z n{z){I + A{z)). 



(1.25) 



The calculation of the principal part of (1.25) at infinity simultaneously determines 
the residues for Y'Y~ l at oo and at 0. The singularity analysis of Y'Y~ X near aj 
for j = 1 , 2 is unchanged from before and one finds 



Y'Y- 1 = -- 
4 





1 



A 



u 



a + v + v 



-u 



+ 



ACi 



XC 2 



z — ai z — a<i 



(1.26) 



where we've written u = u$, v = v$ and v = vq. The analysis of the singularities 
in d a YY -1 proceeds as above and we find (1.15) is unchanged. The non-linear 
equations are deduced as compatibility conditions for (1.26) and (1.15). Once again 
we refer the reader to [14] and [15] for a more detailed look at the resulting equations. 



§2 The tau function 

In this section we will show that the determinant &et(l—XK) which is the focal point 
of the random matrix applications is the r — function for a deformation problem. 
As in the first section there is no serious loss of generality in considering the case of 
just a single interval (ai,a 2 ) and so for simplicity we shall do so. It is also simpler 
to discuss the case in which m(z) = 1, and although it is not essential to do so, we 
will confine our attention in §2 to this special case. 

We begin by reminding the reader of a formula for the logarithmic derivative 
of det(l — \K) which can be found in [15], 



d a logdet(l — XK) = R(ai, ai)da\ — R(d2, a2)dci2- 



(2.0) 



In this formula the terms R(aj,ctj) are the diagonal values for the kernel of the 
resolvent operator, XK(1 — XK) -1 , for XK. In [15] it is shown that 



j-)/ x dq k dp k 
R(a k ,a k ) =Pk— qk- 



dai 



da k 



(2.1) 
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where the notation 

Qk = QM, p k = P(a k ) 7 



is as in 81. We write 



M(z) 



Mu(z) M 12 (z) 
M 21 (z) M 22 (z) 



in (1.16) and find, using (2.1), that, 



R(a k ,a k ) = Tr(M(a k )C k ) + (-1)*"^ 



Tr(dC 2 ) 



(2.2) 



a 2 — a\ 



We write Tr to denote the trace of a matrix. The notation M(cij) is slightly mis- 
leading since M(z) is a polynomial in z with coefficients that are functions of 
a = (ai,a 2 ). To make the dependence on the point a a little clearer we should 
write M(z,a) for M{z). Then M(aj) = M(cij,a). We will use both notations de- 
pending on circumstances. 

Using (2.0) and (2.2) we will make contact with the formula for the log deriva- 
tive of the tau function. In [5] Jimbo, Miwa and Ueno have given a formula for the 
log derivative of the tau function in the general case of irregular singular points. 
In [7] Malgrange has provided a geometrical significance for the tau function in the 
case of regular singular points with simple poles. Instead of simply adopting the 
JMU formula we will introduce the slight modifications in Malgrange's analysis that 
are necessary to deal with the irregular singular point at infinity in (1.14). We hope 
that this will prove to be more instructive. 

First, fix points a° 6 R C C which are distinct for j = 1,2. We suppose that 
the integral operator 1 — K is invertible. Here Kq is the integral operator with 
kernel given by (0.1) and the interval J = (cii,a 2 ). Since 1 — Kq is invertible and 
Kq is compact it is clear that 1 — XKo will remain invertible provided that A G C 
is chosen in a sufficient small neighborhood of 1. Throughout this section we will 
suppose that A is chosen sufficiently close to 1. 

We may rephrase the existence of Y(z,a°) satisfying (1.14) in the following 
terms. There exists a holomorphic connection V° on the trivial bundle C\{a?, a 2 } x 
C 2 with simple poles at a® and a 2 given by, 



V° = d z -Q(z,a°)dz, 



where, 



0(^, a) = M(z,a) + 



ACi(a) AC 2 (a) 



z — 



z — 



and, 



d z = dzd z + d~zd z . 



15 



Note that V° is an integrable connection (it has curvature 0) since 

J^(z,a°)=0 

for z G Cyja^a^}, and locally flat sections for V° are C 2 valued holomorphic 
functions, s(z), satisfying the differential equation 

s'(z) = fl(z,a°)s(z). 

One problem considered by Malgrange in [7] is the problem of deforming the con- 
nection V° in the "a" variables so that the resulting connection, V, has simple poles 
at a,j for j = 1,2 and so that the holonomy of the resulting connection about the 
singular points aj remains fixed as the points cij vary. Actually, to fix the deter- 
mination of holonomy for a connection we must first choose a base point, zq, and 
closed simple paths, 7-,-, about the singular points cij. Parallel translation around 
the curves jj then depends only on the homotopy class of the path 7-,- in C\{ai, CI2}. 
Parallel translation about a closed curve is given by a linear transformation on the 
fiber over the endpoint which our fixed trivialization identifies with C 2 . One might 
consider the problem of deforming the connection V° while keeping the holonomy 
matrices associated with the curves 7^ fixed. However, to make contact with the 
work we did in the first section this is not exactly the problem we will consider. 
Our principal difficulty is that we would like to set zq = 00 G P 1 , but this is not 
a good choice as a base point since the point at infinity is a singular point for the 
connection V. Instead of fixing the holonomy of the connection V we will deform 
V° so that certain monodromy matrices for a distinguished fundamental solution 
Y(z, a) to, 

d z Y(z, a) — Q(z, a)Y(z, a)dz = 0, 

remain fixed as the points a± and a2 vary in a neighborhood of a? and a®. The 
distinguished fundamental solution will arise in the course of fixing the behavior of 
the connection V in a neighborhood of 00. Recall that our starting point in §1 was 
a fundamental solution F(z) to 

d z F - tt(z)Fdz = 0, 

where Q(z) is given by (0.3). We will fix the behavior of V in a neighborhood of 
infinity by requiring that it be gauge equivalent to the connection, 

d z — Vt{z)dz, 

in such a neighborhood, by a gauge transformation which is I + 0(z~ 1 ) near infinity. 
The fundamental solution Y (z, a) will be uniquely determined by the requirement 
that, 

F(z)Y(z,a)- 1 = I + 0(z- 1 ), 



16 



for z in a neighborhood of oo. The monodromy multipliers for Y(z, a) obtained by 
making counterclockwise circuits of a\ and a 2 will play the role of fixed holonomy. 
Analytic continuation of Y(z, a) along a small counterclockwise circuit of a\ will 
transform Y(z,a) into Y(z,a)M, where as in §1, 



Similarly, analytic continuation of Y(z, a) in a small counterclockwise circuit of a 2 
will transform Y(z,a) into Y(z,a)M~ 1 . 

We will next sketch the existence theory for the deformation V by following 
the arguments in Malgrange [7]. In fact, again following [7], it will be convenient 
to do a little more. It is possible to prove the existence of a prolongation of the 
connection V° to a connection in both the z and the a variables in such a fashion 
as to incorporate (1.15) as well as (1.14). 

Let Dj(r) denote the open disk of radius r > about the point a° and write 
Tj(r) for the circle of radius r about the point a°. We suppose that for j = 1,2 real 
numbers pij, p2j, and Tj are chosen so that 



and the numbers p 2 j are small enough so that the closed disks Z?i(p2,i) and 1^2(^2,2) 
do not intersect. We write, 



Choose e > so that e < p\ j for j = 1, 2 and let U denote the open ball of radius 
e about the point (01,02) in C 2 . Let 



We pull back the connection one form Q(z, a°)dz to a connection one form, Vtjd(z — 
cij) on Dj x U\Aj by the translation 



P2,j > Tj > p hj > 



Dj = D j (p 2J ), 
T j = r j( r i) 5 



Aj = {(z,a) E Dj x U\z = dj}. 



(z,a) 



z — a-; + a,. 



One finds 



flj(z, a) = M(z - cij + a°j, a ) + 




z — a* ; + — a 



' 



k 
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where it will be convenient to introduce the notation 

A k (a°) = (-lf-UC^a ). 
One can use this to define a connection 

Vj = d z + d a — Qj(z, a)d(z — cij) 

on the trivial bundle 

(D j x U\Aj) x C 2 . 

Observe that e has been chosen small enough so that the only singularities in Qj(z, a) 
for z G Dj occur when z = aj (where there is a simple pole). Let D denote the 
union of the closed disks Dj(pij) for j = 1,2. Now we extend V° to a connection 
on the trivial bundle 

C\Dx U x C 2 , 

in the following fashion 

Voo = d z + d a - Cl(z, a°)dz. 

Observe that Voo has no singularities on the domain we consider it on. Next we 
will show that the connections Vj and Voo are gauge equivalent to one another on 
the domain SDj x U by a gauge transformation Sj(z, a) which is normalized by the 
condition that Sj(z, a ) = I (the identity matrix). More specifically 

Vj = S^VooSj, (2-3) 

on 5Dj x U. To prove the existence of Sj(z,a) we introduce parallel transport 
relative to the connection Vj. If 7 is a smooth curve in SDj x U then we write 
Pji^f) for parallel transport along the curve 7 relative to the connection Vj. Since 
the bundle we consider is trivial we may think of Pji'j) as a linear map on C 2 . Now 
define 

7j(t, z, a) = (z, (1 — t)a{ + tai, (1 — t)a\ + ta 2 ) for z G Dj, a G U. 

Then 

[0,1] 3t^j(t,z,a), 

is a curve in SDj x U which joins (z,a°) to (z, a). We will denote this curve by 
7j(z, a) and define 

S;\z,a) = Pj( 7 ,M). 
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Observe that we have defined the inverse S 1 " 1 rather than Sj. Because the connec- 
tion Vj is integrable we can replace jj(z, a) by any smooth curve in SDj x U which 
joins (z, a ) to (z, a) but keeps the projection on the first coordinate (i.e., z) fixed. 
The matrix Sj of parallel translation thus satisfies 

d a Sj(z,a) = —Sj{z,a)Vtj{z,a)da. (2-4) 

Now suppose that z G SDj and let Bj(z) be a small ball about z contained in the 
annulus SDj. Suppose that u G Bj(z) and let, 

Pj(u, z\a), 

denote parallel translation with respect to Vj along any smooth curve that joins 
(z, a) to (u, a) whose projection on the first coordinate stays inside Bj(z) and whose 
projection on the second coordinate is fixed at a. Again because the connection is 
integrable and the homotopy class of any such curve is fixed, parallel translation 
does not depend on the choice of such a curve and we have 

d u Pj(u, z\a) = Qj(u, a)Pj(u, z\a)du. (2-5) 

Once again using the fact that parallel translation depends only on the homotopy 
class of the path we can write 

Sj(u, a)Pj(u, z\a) = Pj(u, z\a°)Sj(z, a). 

Differentiating this last relation with respect to u using (2.5) and then setting u = z 
one finds 

d z Sj(z, a) + Sj(z, a)Qj(z, a)dz = Qj(z, a°)Sj(z, a)dz. (2-6) 

But since £lj(z,a°) = Q(z,a°) we can combine (2.4) and (2.6) to get (2.3). 
Now we wish to factor Sj(z, a), 

Sj(z, a) = £00(2, a)~ 1 T,j(z, a), (2.7) 

where Sj(z, a) is a holomorphic, invertible matrix valued function for (z, a) G Dj x 
U, and £00(2:, a) is a holomorphic invertible matrix valued function for z outside 
the union of the disks Dj(pij) and a G U. We further require that £oo(z, a) is 
normalized so that 

EooM = I + 0(z~ 1 ). 

In Malgrange [7] it is shown that by choosing a sufficiently small neighborhood U 
of a this Wiener-Hopf factorization problem can be solved. Suppose now that U 
is sufficiently small so (2.7) is valid. Combining (2.7) and (2.3) we find 

E" 1 ^ = E^VooEoo (2.8) 
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over the set 5Dj x U in the base. On the left hand side of (2.8) is a connection which 
extends to Dj x U\Aj with simple poles at z = aj and on the right hand side the 
connection extends to a connection on C\D x U. We denote by V the connection 
whose existence is thus assured by (2.8). Examining the singularity structure of 
this connection at z = aj and z = oo using the properties of Ej and we find 
that, 

V = d z + d a - M(z, a)dz-J2 ^^-d(z - a k ) , (2.9) 

where M(z, a) is a polynomial in z with the same leading terms as M(z, a ) and 
hence the same leading terms as Q(z) near infinity, and A^{a) agrees with Ak(a°) 
at a = a . To make contact with the developments in §1 we introduce Cfc(a) by 

A k (a) = (-l^XCkia). 

The distinguished fundamental solution for V whose monodromy is to remain fixed 
as a function of a is given by, 

Y(z,a) = ^ 1 (z,a)Y(z,a°), 

for (z, a) G C\DxU. It is clear that this has exactly the same monodromy multiplier 
as Y(z, a ), as z makes a circuit of either of the circles Tj. 

Now we turn to the definition of the r — function. Let H denote the Hilbert 
space, 

H = H*(T 1 )@H*(T 2 % 

where denotes the Sobolev space of order ^. Let H + denote the subspace of 
H which consists of boundary values of functions analytic in the open disks Dj(rj) 
for j = 1,2 (recall that the circle 1^ is the boundary of the disk Dj(rj)). Let H- 
denote the boundary values in H of functions analytic in 

which vanish at z = oo. Then H = H + ©i7_ and we denote the projections relative 
to this decomposition by P±. If S is a smooth matrix valued function defined on 
the circles Fj for j = 1, 2 it defines a map 

acting by multiplication by S in the obvious way. The Toeplitz operator, T$, asso- 
ciated with S is defined by 

T s = P+SP+. 
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We are now prepared to define the r — function. Let S denote the matrix valued 
function which is equal to Sj on Tj for j = 1,2. Then, following Malgrange, we 
define 



The inverse T s 1 is computable in terms of the factorization (2.7) by the standard 
Wiener-Hopf technique and (see lemma 6.7 in [7]) one finds 



In [10] the r — function defined by (2.10) is shown to have the following significance. 
The family of subspaces 



determine boundary conditions for a Cauchy-Riemann operator localized in the 
exterior of D. This Cauchy-Riemann operator has a domain which incorporates 
functions which have suitable branching at the points aj. The subspaces W(a) 
belong to a Grassmannian over which there sits a holomorphic line bundle det*. 
This line bundle has a connection whose curvature vanishes over the family of 
subspaces W(a). The r — function defined by (2.10) is the ratio of the canonical 
section of det* to a section of det* that is flat with respect to this connection. It is 
also explained in [10] that this makes it possible to regard r as the determinant of 
a singular Cauchy-Riemann operator. The reader is refered to [10] for details. 
From (2.8) and (2.9) we see that since 



dlogr = Tr (Tg 1 Td a s - T s -i daS ) ■ 



(2.10) 




j 



(2.11) 



W(a) = Y(-,a°)- 1 S(-,a)H+ 




we have 




(2.12) 



Using (2.8) and (2.9) and 



V, 



one sees that 



Ej 1 (z,a)d z Y, j (z,a) 




(2.13) 



where 
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Since T,j(z,a) is an invertible holomorphic function on Dj x U it follows that the 
pole type singularities in the two terms on the right hand side of (2.13) must cancel. 
One may easily confirm that this implies 

A j (a) = a- 1 A j (a°)a j , (2.14) 

where 

dj := T, j (a j ,a). 

We write 

Zj{z, a) = a,(I + f3j{z - a,)) + O ((z - a,) 2 ) 

for the first terms in the Taylor expansion of S 7 - near z = ctj. Using (2.13) and 
(2.14) one finds that 

E-^E,- = B 3 {a) - a' 1 (^(a ) - [A^aP), (3 3 \) a 3 +0((z- a,)) , (2.15) 

where 



B j (a) = M(a j ,a) + J2- 



A k (a) 



Substituting (2.15) and (2.12) into the formula (2.11) one finds 

d a logr = J^Tr (Bji^Ajia) - Bj{a )Aj (a )) daj, (2.16) 

3 

where we made use the centrality of the trace, (2.14) again, and the fact that 

Tr([A j (a%0 j ]A j (a°))=O. 
If we use Aj(a) = (-l) j - 1 XC j (a) to compare (2.16) with (2.0) and (2.2) we find 

d a logdet(l - \K) =^2Tr(B j (a)A j (a))da j , 

3 

from which it follows that the r— function differs from det(l— \K) by the exponential 
of a linear function of a. 
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